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c^ . Abstract 

In this paper, the hnear free flexural vibrations of functionally graded mate- 
rial plates with a through center crack is studied using an 8-noded shear flex- 
^ ' ible element. The material properties are assumed to be temperature depen- 

^ . dent and graded in the thickness direction. The effective material properties 

Q ! are estimated using the Mori-Tanaka homogenization scheme. The formula- 

'vh I tion is developed based on first-order shear deformation theory. The shear 

t^^ ' correction factors are evaluated employing the energy equivalence principle. 

■ The variation of the plates natural frequency is studied considering various 

parameters such as the crack length, plate aspect ratio, skew angle, temper- 
ature, thickness and boundary conditions. The results obtained here reveal 
'k> I that the natural frequency of the plate decreases with increase in temperature 

^ ' gradient, crack length and gradient index. 

C^ ! 
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1. Introduction 

A functionally graded material (FGM) is a new class of material whose 
properties are characterized by the volume fraction of its constituent mate- 
rials. The concept of characterization of material properties is not new, but 
the unique feature of FGMs is that these materials are made of a mixture 
of ceramics and metals that are characterized by the smooth and continuous 
variation in properties from one surface to another [l|, Ij, y] . For structural 
integrity, FGMs are preferred over fiber-matrix composites that may result 
in debonding due to the mismatch in the mechanical properties across the 
interface of two discrete materials bonded together. With the increased use 
of these materials for structural components in many engineering applica- 
tions, it is necessary to understand the dynamic characteristics of function- 
ally graded plates. 

Many researchers have recently attempted to study the bending behav- 
ior of FGM plates on three-dimensional elasticity solutions j^, |5|, Ig, lZJ . All 
these works are limited to simply supported plates under sinusoidal trans- 
verse mechanical or thermal loading. Reddy and Cheng [^ and Vel and 
Batra [5| have accounted for the variation of material properties through 
the thickness according to a power-law distribution and the locally effective 
material properties were obtained in terms of the volume fractions of the con- 
stituents through the Mori-Tanaka homogenization scheme. Kashtalayan J6| 
derived the elasticity solutions making use of the Plevako general solution 
of the equilibrium equations for inhomogeneous isotropic media, whereas. 
Pan [Tl] studied the laminated functionally graded simply supported rect- 
angular plates under sinusoidal surface load, extending the Pagano's solu- 
tions which may not be valid for finding the solutions of such plate problems 
with continuous inhomogeneity. Elishakoff and Gentilini [8|] investigated the 
three-dimensional static analysis of clamped functionally graded plates under 
uniformly distributed load applying the Ritz energy method. 

Application of 3D analysis, in general, is quite cumbersome while deal- 
ing with complex loading and boundary conditions. Hence, the analysis of 
isotropic, composite and FGM plates is carried out numerically as well as an- 
alytically using plate theories assuming plane stress conditions. Such approx- 
imation can predict global displacement and bending moments with sufficient 
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accuracy [9|, |lOJ. Few analytical and finite element studies on the bending 
analysis of FGM plates are recently available in the literature using plate the- 
ories. Qian jll| and Matsunaga [12| examined the bending of thick square 
FGM plates considering the higher-order shear deformation theory, whereas, 
Zenkour 13|, |lj] dealt with 2D trigonometric functions based shear deforma- 
tion theory. The nonlinear thermo-mechanical response of FGM plates was 
examined by Praveen and Reddy 15| and Reddy 16|] considering the higher- 
order structural theory. Carrera ITJ have obtained closed form and finite ele- 
ment solutions for the static analysis of functionally graded plates subjected 



to transverse mechanical loads. The unified formulation employed in [17 
permits a large variety of plate models with variable kinematic assumptions 
covering first-order as well as higher-order theories. However, higher-order 
models that involve additional displacement fields may be based on either an 
equivalent single layer theory or discrete layer approach. Furthermore, they 
are computationally expensive in the sense that the number of unknowns 
to be solved is high compared to that of the first-order shear deformation 
formulation. 

In literature, there has been many studies on dynamic characteristics of 



FGM plates [15|, ll8|, ll9|, |20|, l2l|, |22|, |23|, Mi and shells |25|, |26|, |27|, l28|, |29 



The dynamic characteristics of a cracked structural element is especially im- 
portant because a crack in a vibrating structure results in stiffness decrease, 
stress concentration, anisotropy and local flexibility, which are functions of 
the crack location and size. Moreover the crack will open and close depend- 
ing on the vibration amplitude. The vibration of cracked isotropic plates 
was studied as early as 1969 by Lynn and Kumbasa r |30l| who used a Green's 
function approach. Later, in 1972, Stahl and Keer |3l| studied the vibration 
of cracked rectangular plates using elasticity methods. The other numer- 
ical methods that are used to study the dynamic response and instability 
of plates with cracks or local defects are: (1) Finite fourier series trans- 
form 32|; (2) Rayleigh-Ritz Method (33|; (3) harmonic balance method 34 1; 
(4) finite element method [35|,|36|; (5) extended finite element method |37] : 



(6) Smoothed finite element methods 
methods 



39, 40, 41 



43, 11, 44, 45 



42 1 and (7) Meshfree 
Recently, Yang and Chen 46[ and Kitipornchai 
et al., 47| studied the dynamic characteristics of FGM beams with an edge 
crack. In case of FGM, the dynamic characterisitcs also depend on the gra- 
dient index compared to the isotropic case. However, to the author's knowl- 
edge, studies of cracked FGM plates are scarce in the literature that are of 
practical importance to the designers. 



The main focus of this paper is to compute the hnear free flexural vi- 
brations of FGM plates with a through center crack using the finite element 
method. Here, an eight-noded shear flexible quadrilateral plate element based 



on field consistency approach |48l . |49| is used to study the dynamic charac- 



teristics of FGM plates subjected to thermo-mechanical loadings. The shear 
correction factor is calculated from the energy equivalence principle. The 
temperature field is assumed to be constant in the plate and varied only in 
the thickness direction. The material is assumed to be temperature depen- 
dent and graded in the thickness direction according to a power law distri- 
bution in terms of the volume fractions of the constituents. The effective 
material properties are estimated from the volume fractions and the ma- 
terial properties of the constituents using the Mori-Tanaka homogenization 



method [50, |511 ■ The formulation developed herein is validated considering 
different problems for which the solutions are available in the literature. De- 
tailed parametric studies are carried out to understand the free vibration 
characteristics of cracked FGM plates. 

The paper is organized as follows. In Section [21 a brief introduction to 
FGM is given followed by the plate formulation. The treatment of boundary 
conditions for skew plates and eight-noded quadrilateral plate element are 
discussed in Section [31 A detailed numerical study is presented in Section [H 
followed by conclusion in the last section. 

2. Theoretical development and formulation 

2.1. Functionally Graded Material 

A functionally graded material (FGM) rectangular plate (length a, width 
b and thickness h), made by mixing two distinct material phases: a metal 
and ceramic is considered with coordinates x, y along the in-plane directions 
and z along the thickness direction (see Figure dll)). The material on the 
top surface {z = h/2) of the plate is ceramic and is graded to metal at 
the bottom surface of the plate {z = —h/2) by a power law distribution. 
The homogenized material properties are computed using the Mori-Tanaka 



Scheme 50, 51 



2.1.1. Estimation of mechanical and thermal properties 

Based on the Mori-Tanaka homogenization method, the effective bulk 
modulus K and shear modulus G of the FGM are evaluated as 50|, [5l|, l52|, [ll 
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Here, Vi {i = c, m) is the volume fraction of the phase material. The sub- 
scripts c and m refer to the ceramic and metal phases, respectively. The 
volume fractions of the ceramic and metal phases are related by Vc + Kn = 1, 
and Vc is expressed as 
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where n in Equation (^ is the volume fraction exponent, also referred to as 
the gradient index. The effective Young's modulus E and Poisson's ratio u 
can be computed from the following expressions: 
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The effective heat conductivity coefficient k and coefficient of thermal ex 
pansion a is given by |53|, |54 



K — 


K„iYi 


Kc- 


^m 


a — 


am 



ar 



OLr. 



Vr 



1 + (1 - K 






Kr, 



K„ 



(5) 



The effective mass density p is given by the rule of mixtures as [20 



P = PcVc + pmVm 



(6) 



2.1.2. Temperature distribution through the thickness 

Th e material properties P that are temperature dependent can be written 
as 
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P = PoiP-iT-^ + 1 + PiT + PsT' + PsT^), (7) 

where Po, P-i, Pi, P2, P3 are the coefficients of temperature T and are unique 
to each constituent material phase. The temperature variation is assumed to 
occur in the thickness direction only and the temperature field is considered 
to be constant in the xy-plane. In such a case, the temperature distribution 
along the thickness can be obtained by solving a steady state heat transfer 
equation 



d 

dz 



' .dT 



0, T = Tc at z = h/2- T = T^at z = -h/2 



The solution of Equation (|8]) is obtained by means of a polynomial series [56 
as 
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where K^m = i^c~ i^m and Tc, T^ denote the temperature of the ceramic and 
metal phases, respectively. 




(a) 




(b) 

Figure 1: (a) coordinate systeni of a rectangular FGM plate, (b) Coordinate system of a 
skew plate 

2. 2. Plate formulation 

Using the Mindlin formulation, the displacements u, v,w ai a point (x, y, z) 
in the plate (see Figure ([1])) from the medium surface are expressed as func- 
tions of the mid-plane displacements Mq, Wq, Wq and indepedent rotations 6^, Oy 
of the normal in yz and xz planes, respectively, as 



u{x,y,z,t) 
v{x,y,z,t) 
w{x,y,z,t) 



Uo{x,y,t) + z6^{x,y,t) 
Vo{x,y,t) + z9y{x,y,t) 
Wo{x,y,t) 



(12) 



where t is the time. The strains in terms of mid-plane deformation can be 



written as 
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(13) 



The midplane strains £p, bending strain 8^, shear strain Eg in Equation ( 1T5]) 
are written as 
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(14) 



where the subscript 'comma' represents the partial derivative with respect to 
the spatial coordinate succeeding it. The membrane stress resultants N and 
the bending stress resultants M can be related to the membrane strains, Sp 
and bending strains Sb through the following constitutive relations 
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where the matrices a = Ajj,B = Bij and D = Dij] {i,j = 1,2,6) are the 
extensional, bending-extensional coupling and bending stiffness coefficients 
and are defined as 
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The thermal stress resultant, N^ and the moment resultant M'^ are 



(16) 



XX 



yy ( / ^v 

xy ' 



h/2 



^ = < K }= Q.Mz,T){ 1 } AT{z)dz 



m: 



XX 

T 



h/2 
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where the thermal coefficient of expansion a{z,T) is given by Equation (^ 
and AT{z) = T{z)—To is the temperature rise from the reference temperature 
To at which there are no thermal strains. 

Similarly, the transverse shear force Q = {Qxz,Qyz} is related to the 
transverse shear strains Eg through the following equation 

Qij = EijEs (18) 

where Eij = /_)j/2 Q'^i'^i dz; {i,j = 4,5) is the transverse shear stiffness 
coefficient, Vi,Vj is the transverse shear coefficient for non- uniform shear 
strain distribution through the plate thickness. The stiffness coefficients Q^^j 
are defined as 



7^ _7T _Eiz,T) - _ uEiz,T ) 



Qll — Q22 ~ 1 ,,2 ' ^12 ~ . 2 ' ^16 ~ ^^26 ~ U 



where the modulus of elasticity E{z,T) and Poisson's ratio u are given by 
Equation (jlj). The strain energy function U is given by 



U{S) = - / {e'pAsp + e'pBsb + sjBsp + s'J^Sb + e^Es 
^ Jn 
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(20) 
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where 6 = {u,v,w,9x,9y} is the vector of the degree of freedom associated 
to the displacement field in a finite element discretization. Following the 
procedure given in [57|, the strain energy function U given in Equation fl20l) 
can be rewritten as 

9 



U{S) = ^<5^K5 
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where K is the hnear stiffness matrix. The kinetic energy of the plate is 
given by 
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where p = /_^ ^ p{z) dz, I = J_^'j^% z^p{z) dz and p{z) is the mass density 
that varies through the thickness of the plate given by Equation ([6]). The 
plate is subjected to a temperature field and this in turn results in in-plane 
stress resultants (A^^^, Nyy, N^^). The external work due to the in-plane stress 
resultants (A^^^, Nyy, A^j;.))) developed in the plate under the thermal load is 



V{S) 




Nl'tf^] +Ar:M^) +2Ar:M 
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(23) 



Substituting Equation ( 12T|) - (123!) in Lagrange's equations of motion, the 
following governing equation is obtained 



M5 + (K + Kg)5 = (24) 

where M is the consistent mass matrix. After substituting the characteris- 
tic of the time function |49| 6 = —u'^S, the following algebraic equation is 
obtained 



((K + Kg) - uj'^M) 6 = (25) 

where K and Kg are the stiffness matrix and geometric stiffness matrix 
due to thermal loads, respectively, u is the natural frequency. The plate is 
subjected to a temperature field. So, the first step in the solution process is 
to compute the in-plane stress resultants due to the temperature field. These 
will then be used to compute the stiffness matrix of the system and then the 
frequencies are computed for the system. 
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3. Element Description and Treatment of boundary conditions 

3.1. Eight noded shear flexible element 

The plate element employed here is a Cq continuous shear flexible element 
with flve degrees of freedom {uo,Vo,Wo,0x,9y) at eight nodes in an 8-noded 
quadrilateral (QUAD-8) element. If the interpolation functions for QUAD-8 
are used directly to interpolate the flve variables {uo, Vo, Wo, 0^, By) in deriving 
the shear strains and membrane strains, the element will lock and show os- 
cillations in the shear and membrane stresses. The fleld consistency requires 
that the transverse shear strains and membrane strains must be interpolated 
in a consistent manner. Thus, the Ox and 9y terms in the expressions for 
shear strain Ss have to be consistent with the derivative of the fleld func- 
tions, Wo^x and Wo^y This is achieved by using fleld redistributed substitute 
shape functions to interpolate those speciflc terms, which must be consistent 



as described in [48|, |49|. This element is free from locking syndrome and 



has good convergence properties. For complete description of the element. 



interested readers are referred to the literature |48l . |49| , where the element 
behavior is discussed in great detail. 

3.2. Skew boundary transformation 

For skew plates supported on two adjacent edges, the edges of the bound- 
ary elements may not be parallel to the global axes {x,y,z). In order to 
specify the boundary conditions at skew edges, it is necessary to use edge 
displacements {u'^, v'^, w'^) , etc., in local coordinates {x',y',z') (see Figure 
([T])).The element matrices corresponding to the skew edges are transformed 
from global axes to local axes on which the boundary conditions can be con- 
veniently specifled. The relation between the global and local degrees of 
freedom of a particular node can be obtained through simple transformation 



rules [58| expressed as 



d = Lgd' (26) 

where di, d'^ are the generalized displacement vectors in the global and local 
coordinate system, repectively, of a node i and are deflned as 

d = [Uo Vo Wo Ox Oy ] 

d' = [u'o < < 9'^ e'^Y (27) 
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The nodal transformation matrix for a node i, on the skew boundary is given 
by 



cos IJJ 


sin ijj 








- sin ijj 


cos ijj 











1 














cosip 


sinip 











— sin ip cos tp 



(28) 



where ip is the angle of the plate. The transformation matrix for the nodes 
that do not lie on the skew edge is an identity matrix. Thus, for the complete 
element, the element transformation matrix is written as 



diag( L„ L, L„ L„ L„ L, 



L. > 



(29) 



4. Numerical Examples 

In this section, we present the natural frequencies of cracked FGM plates. 
Figure (|2]) shows the geometry and boundary condition of the plate. The 
crack is assumed to be at the center of the plate. For a rectangular plate, 
the crack is assumed to be parallel to the sides of length a. The effect of 
the gradient index n, temperature gradient (Tc — T^), crack length c, aspect 
ratio a/b, thickness of the plate h, skew angle of the plate ip and boundary 
conditions on the natural frequency are numerically studied. The assumed 
values for the parameters are listed in Table [TJ 

In all cases, we present the non-dimensionalized free flexural frequency 
defined as 



Vt = ua\l^ (30) 

where u is the natural frequency, D = ^20-1^2) ^^ ^^^ bending rigidity of a 
homogeneous ceramic plate, pc, E^ are the mass density and Young's modulus 
of the ceramic, respectively. Figure (|3]) shows the variation of the volume 
fractions of ceramic and metal, respectively, in the thickness direction z for 
the FGM plate. The top surface is ceramic rich and the bottom surface is 
metal rich. The FGM plate considered here consists of silicon nitride (SI3N4) 
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(a) Simply supported (b) clamped 

Figure 2: Plate with a center crack 



Table 1: Details of parameters used in the numerical study 



Parameter 



Assumed values 



Gradient index, n 

Ceramic Temperature (K), T^ 

Metal Temperature (K), T^ 

Crack Length c/a 

Aspect ratio a/h 

Thickness of the plate a/h 

Skew angle of the plate if) (in deg.) 

Boundary conditions 



0, 0.5, 1,2, 5, 10 

300, 400, 600 

300, 300, 300 

0.2, 0.4, 0.6 

1,2 

10, 20 

0°, 15°, 30°, 45° 

Simply supported and Clamped 
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gradient index, n 
n=0.5 — ^ 

n=1 
n=2 
n=5 
n=10 




Figure 3: Through thickness variation of volume fraction 
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Table 2: Temperature dependent coefficient for material SI3N4/SUS3O4, Ref. |55l.l21| 



Material Property Po P_i Pi P2 



SI3N4 



SUS304 



E(Pa) 348.43e9 
a (1/K) 5.8723e-*^ 



E(Pa) 201. 04e^ 
a (1/K) 12.330e-*^ 



0.0 


-3.070e-^ 


2.160e"^ 


-8.946e-i^ 


0.0 


9.095e-4 


0.0 


0.0 


0.0 


3.079e-^ 


-6.534e-^ 


0.0 


0.0 


8.086e-4 


0.0 


0.0 



and stainless steel (SUS304). The material is considered to be temperature 
dependent and the t emp erature coefficients corresponding to SI3N4/SUS304 



are listed in Table |2] [2l|, |55| . The mass density (p) and thermal conductivity 
(K) are: Pc=2370 kg/m^, Kc=9.19 W/mK for SI3N4 and Pm = 8166 kg/m^, 
Km = 12.04 W/mK for SUS304. Poisson's ratio u is assumed to be constant 



and taken as 0.28 for the current study |21l . |22j. Here, the modified shear 



correction factor obtained based on energy equivalence principle as outlined 



in [59| is used. The boundary conditions for simply supported and clamped 

cases are (see Figure ([2])): 

Simply supported boundary condition: 



Uo = Wo = Oy = on X = 0,a 

Vo = Wo = Or, = ony = 0,b (31) 

Clamped boundary condition: 

Uo = Wo = Oy = Vo = Ox = on X = 0, a 

Uo = Wo = Oy = Vo = 0^ = ony = 0,b (32) 

Before proceeding with the detailed study on the effect of different parame- 
ters on the natural frequency, the formulation developed herein is validated 
against available results pertaining to the linear frequencies for functionally 



graded plate based on higher-order shear deformation plate theory |44j . Ill 



and to the linear frequencies for cracked isotropic plates [60|, 161|- The com- 
puted frequencies for a square simply supported plate: (a) with gradient 
index n =1 and a/h = 10 is given in Table [31 (b) with a center crack and 
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thickness a/h = 1000 are given in Table HI It can be seen that the numerical 
results from the present formulation are found to be in good agreement with 
the existing solutions. 

It is seen that with increase in crack length c/a, the natural frequency 
decreases. The decrease in frequency with increasing crack length c/a is due 
to reduction in the stiffness of the plate. The linear frequencies for uncracked 
FGM plates in thermal environment is shown in Table and the results are 
compared with the results available in the literature [62[. Here the calculated 
non-dimensional frequency is defined as: fi = uj{a^ /h)^J pm{^ — y'^)/ Em- 
Based on the progressive mesh refinement, an 8x8 mesh is found to be ade- 
quate to model the full plate for the uncracked FGM plates. The mesh size 
used for cracked FGM plates and number of degrees of freedom are given in 
Table H 

Table 3: Non-dimensionalized natural frequency VI — ujh^J pm/ Em for a simply supported 
square plate for_? = 10, Gradient index, n—1. The material properties are assumed to be 



the same as in [4J, [11| . The present formulation based on First order Shear Deformation 
plate Theory is compared with Third order Shear Deformation Theory (TSDT) [4^] and 
Higher Or der Shear and Normal Deformable Plate Theory (HOSNDPT) ^. 

Mode TSDT [4^ HOSDNPT [11] FSDT (Present) 

1 0.0592 0.0584 0.0609 

2 0.1428 0.1410 0.1422 

3 0.1428 0.1410 0.1422 



Next, the linear flexural vibration behavior of FGM is numerically studied 
with and without thermal environment. For the uniform temperature case, 
the material properties are evaluated at T^. = T^, = T = 300K. The variation 
of natural frequency with aspect ratio, thickness, gradient index and crack 
length is shown in Table [7] and graphically in Figures dl]) - dZ]) . Figure (jl]) 
and Figure (|5]) shows the variation of mode 1 and mode 2 frequency of a 
simply supported square plate with gradient index for different crack length 
c/a and thickness h. Figure (j6]) and Figure ([7]) shows the variation of mode 1 
and mode 2 frequency of a simply supported rectangular plate with gradient 
index for different crack length c/a and thickness h. The natural frequency 
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Table 4: Non-dimensionalized natural frequency for a simply supported square plate with 
a center crack for f = 1000. 



c/a 



Mode 1 



Mode 2 



Ref [60] Ref [6JJ Present Ref [60] Ref [61] Present 



0.0 



19.74 



19.74 19.752 



49.35 49.38 49.410 



0.2 



19.38 19.40 19.399 



49.16 49.85 49.300 



0.4 



18.44 18.50 18.484 



46.44 47.27 47.01^ 



0.6 



17.33 17.37 17.439 



37.75 38.92 38.928 



Table 5: Comparison of non-dimensional linear frequencies of simply supported FGM plate 
(a/6 = 1, a/h — 8) in thermal environment. 



Temperature 


gradient 
n 


Mode 1 


Mode 2 


T T 


Ref. [62 


Present 


Ref. [62 


Present 




0.0 


12.397 


12.311 


29.083 


29.016 


T™=300, 

re=4oo 


0.5 


8.615 


8.276 


20.215 


19.772 




1.0 


7.474 


7.302 


17.607 


17.369 




2.0 


6.693 


6.572 


15.762 


15.599 




0.0 


11.984 


11.888 


28.504 


28.421 


7^=300, 
Tc=600 


0.5 


8.269 


7.943 


19.784 


19.327 




1.0 


7.171 


6.989 


17.213 


16.959 




2.0 


6.398 


6.269 


15.384 


15.207 
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Table 6: Number of nodes and total number of degrees of freedom for FGM plates with 
and without crack used in the present analysis. 



c/a 


Number 


Total degrees 




of nodes 


of freedom 


0.0 


225 


513 


0.2 


184 


840 


0.4 


264 


1216 


0.6 


356 


1604 



increases with decrease in the plate thickness and this behavior is maintained, 
regardless of the boundary conditions. It is seen from Table [7] and Figures 
(jlj) - ([7]) that with increase in crack length and gradient index, the linear 
frequency of the plate decreases and with decrease in plate thickness, the 
natural frequency increases. The decrease in frequency with increasing crack 
length cja is due to the reduction in the stiffness of the plate and the decrease 
in frequency with increasing gradient index n occurs due to increase in the 
metallic volume fraction. In both the cases, the reduction in frequency is 
due to the stiffness degradation. It can also be observed that the decrease in 
frequency is not proportional to the increase in crack length. The frequency 
decreases rapidly with increase in crack length. This observation is true, 
irrespective of thick or thin, and square or rectangular case considered in the 
present study. 
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Figure 4: Frequency (Mode 1) as a function of gradient index n and crack length c/a for 
a simply supported square FGM plate in ambient temperature {Tc ~ 300K, Tm ~ 300K). 



19 



T3 

o 



D" 



T3 
O 



45 



40 



35 



30 



25 



20 



15 



10 



\ 






Crack Length, c/a 


c/a=0.2 - -X- -- _ 
c/a=0.4 >^ 
c/a=0.6 B 


' A 








■■■ ^, \_ 








V, 


^^^^^^^^4^ 




- 




"*'""""■■■------.. 


^^'^'^'^'^^-^*— 




''""""'■---><- 




- 




□ 


-^ 



4 6 

Gradient index, n 



10 



(a) a/h=lO 



Crack Length, c/a 
c/a=0.0 — ^ 

c/a=0.2 --X- 

c/a=0.4 » 

c/a=0.6 




4 6 

Gradient index, n 



10 



(b) a/h^20 

Figure 5: Frequency (Mode 2) as a function of gradient index n and crack length c/a for 
a simply supported square FGM plate in ambient temperature {Tc ~ 300K, Tm ~ 300K). 
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Figure 6: Frequency (Mode 1) as a function of gradient index n and crack length c/a for 
a simply supported rectangular FGM plate in ambient temperature (Tc = 300K, T™ = 
300K). 
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(b) a/h=20 

Figure 7: Frequency (Mode 2) as a function of gradient index n and crack length c/a for 
a simply supported rectangular FGM plate in ambient temperature (Tc = 300K, T™ = 
300K). 
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(b) Mode 2 

Figure 8: Frequency as a function of gradient index n and different temperature gradient 
for a simply supported rectangular FGM plate. The results are shown for crack length 
c/a = 0.2. 
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Table 7: Effect of aspect ratio, thickness, crack length and gradient index on the natural frequency of simply supported FGM 
plates in ambient temperature {Tc ~ 300K, Tm = 300K). 





a/b a/h 


n 




Mode 1 






Mode 2 








c/ 


'a 






c/ 


'a 






0.0 


0.2 


0.4 


0.6 


0.0 


0.2 


0.4 


0.6 






0.0 


18.357 


17.815 


16.805 


15.829 


43.923 


43.316 


39.765 


31.671 






0.5 


12.584 


12.214 


11.522 


10.853 


30.032 


29.622 


27.218 


21.697 




10 


1.0 


11.069 


10.742 


10.132 


9.544 


26.371 


26.008 


23.888 


19.036 




2.0 


9.954 


9.659 


9.111 


8.582 


23.703 


23.372 


21.447 


17.077 






5.0 


9.026 


8.759 


8.262 


7.782 


21.527 


21.222 


19.455 


15.479 




1 


10.0 


8.588 


8.334 


7.861 


7.405 


20.511 


20.221 


18.538 


14.749 




0.0 


18.829 


18.340 


17.366 


16.377 


46.599 


46.178 


43.337 


35.265 






0.5 


12.908 


12.573 


11.906 


11.228 


31.862 


31.575 


29.643 


24.133 


to 


20 


1.0 


11.353 


11.058 


10.470 


9.873 


27.983 


27.730 


26.028 


21.184 




2.0 


10.212 


9.946 


9.416 


8.879 


25.165 


24.936 


23.396 


19.033 






5.0 


9.263 


9.021 


8.541 


8.054 


22.868 


22.658 


21.252 


17.281 






10.0 


8.813 


8.584 


8.127 


7.664 


21.786 


21.586 


20.247 


16.465 






0.0 


43.821 


39.836 


33.381 


28.147 


67.409 


67.127 


65.149 


48.012 






0.5 


30.052 


27.315 


22.877 


19.282 


46.123 


45.930 


44.576 


32.901 




10 


1.0 


26.433 


24.023 


20.116 


16.953 


40.510 


40.341 


39.153 


28.867 




2.0 


23.761 


21.595 


18.085 


15.245 


36.397 


36.246 


35.181 


25.888 






5.0 


21.532 


19.572 


16.399 


13.830 


33.026 


32.889 


31.924 


23.450 




2 


10.0 


20.485 


18.622 


15.606 


13.163 


31.460 


31.329 


30.409 


22.342 




0.0 


46.468 


42.594 


35.886 


30.178 


73.591 


73.283 


71.260 


56.844 






0.5 


31.864 


29.205 


24.596 


20.676 


50.361 


50.151 


48.766 


38.920 




20 


1.0 


28.031 


25.684 


21.621 


18.172 


44.251 


44.067 


42.846 


34.154 




2.0 


25.211 


23.095 


19.439 


16.338 


39.792 


39.626 


38.526 


30.667 






5.0 


22.860 


20.943 


17.634 


14.826 


36.131 


35.980 


34.982 


27.826 






10.0 


21.747 


19.927 


16.784 


14.114 


34.408 


34.264 


33.315 


26.513 



Table 8: Effect of aspect ratio, thickness, crack length and gradient index on the natural frequency of simply supported FGM 
plates with a temperature gradient {Tc = 400K, Tm = 300K) 





a/b a/h 


n 




Mode 1 






Mode 2 








c/ 


'a 






c/ 


'a 






0.0 


0.2 


0.4 


0.6 


0.0 


0.2 


0.4 


0.6 






0.0 


17.962 


17.433 


16.465 


15.537 


43.382 


42.803 


39.331 


31.338 






0.5 


12.281 


11.920 


11.261 


10.630 


29.640 


29.251 


26.907 


21.458 




10 


1.0 


10.786 


10.468 


9.889 


9.337 


26.014 


25.672 


23.606 


18.820 




2.0 


9.862 


9.397 


8.877 


8.383 


23.366 


23.055 


21.183 


16.876 






5.0 


8.753 


8.495 


8.028 


7.584 


21.196 


20.912 


19.199 


15.285 




1 


10.0 


8.309 


8.065 


7.623 


7.204 


20.177 


19.909 


18.282 


14.556 




0.0 


17.533 


17.095 


16.275 


15.463 


45.169 


44.850 


42.272 


34.474 






0.5 


11.880 


11.585 


11.041 


10.505 


30.756 


30.552 


28.827 


23.529 


to 


20 


1.0 


10.379 


10.122 


9.651 


9.189 


26.947 


26.772 


25.267 


20.622 




2.0 


9.255 


9.026 


8.612 


8.209 


24.158 


24.006 


22.660 


18.491 






5.0 


8.276 


8.073 


7.713 


7.365 


21.843 


21.714 


20.507 


16.735 






10.0 


7.785 


7.596 


7.266 


6.949 


20.730 


20.614 


19.483 


15.905 






0.0 


43.280 


39.345 


32.999 


27.855 


66.733 


66.461 


64.522 


47.643 






0.5 


29.651 


26.951 


22.596 


19.069 


45.642 


45.456 


44.132 


32.647 




10 


1.0 


26.067 


23.690 


19.860 


16.760 


40.077 


39.915 


38.754 


28.644 




2.0 


23.416 


21.281 


17.844 


15.064 


35.995 


35.849 


34.810 


25.685 






5.0 


21.196 


19.266 


16.166 


13.657 


32.638 


32.507 


35.168 


23.262 




2 


10.0 


20.149 


18.317 


15.374 


12.992 


31.073 


30.949 


30.057 


22.160 




0.0 


45.036 


41.302 


34.921 


29.480 


72.011 


71.734 


69.833 


56.154 






0.5 


30.755 


28.204 


23.847 


20.137 


49.155 


48.968 


47.677 


38.412 




20 


1.0 


26.992 


24.746 


20.921 


17.667 


43.129 


42.965 


41.833 


33.691 




2.0 


24.202 


22.184 


18.760 


15.851 


38.710 


38.565 


37.552 


30.230 






5.0 


21.835 


20.018 


16.948 


14.337 


35.044 


34.914 


34.006 


27.400 






10.0 


20.691 


18.974 


16.079 


13.614 


33.295 


33.175 


32.320 


26.087 



Table 9: Effect of aspect ratio, thickness, crack length and gradient index on the natural frequency of simply supported FGM 
plates with a temperature gradient {Tc = 600K, Tm = 300K) 





a/b a/h 


n 




Mode 1 






Mode 2 








c/ 


'a 






c/ 


'a 






0.0 


0.2 


0.4 


0.6 


0.0 


0.2 


0.4 


0.6 






0.0 


17.105 


16.605 


15.730 


14.910 


42.272 


41.754 


38.451 


30.666 






0.5 


11.611 


11.273 


10.687 


10.142 


28.808 


28.467 


26.252 


20.958 




10 


1.0 


10.155 


9.859 


9.350 


8.878 


25.241 


24.944 


23.000 


18.357 




2.0 


9.064 


8.799 


8.349 


7.935 


22.617 


22.350 


20.596 


16.429 






5.0 


8.115 


7.788 


7.483 


7.123 


20.425 


20.187 


18.598 


14.828 




1 


10.0 


7.642 


7.420 


7.054 


6.723 


19.371 


19.150 


17.654 


14.079 




0.0 


14.325 


14.017 


13.611 


13.269 


41.964 


41.849 


39.931 


32.750 






0.5 


9.269 


9.082 


8.881 


8.734 


28.239 


28.165 


27.002 


22.190 


to 


20 


1.0 


7.866 


7.712 


7.573 


7.492 


24.564 


24.501 


23.545 


19.360 




2.0 


6.725 


6.602 


6.531 


6.515 


21.812 


21.759 


20.971 


17.254 






5.0 


5.543 


5.458 


5.483 


5.565 


19.401 


19.358 


18.755 


15.455 






10.0 


4.818 


4.761 


4.864 


5.023 


18.161 


18.125 


17.645 


14.564 






0.0 


42.172 


38.336 


32.221 


27.270 


65.381 


65.129 


63.276 


46.934 






0.5 


28.805 


26.181 


22.005 


18.626 


44.640 


44.470 


43.210 


32.132 




10 


1.0 


25.277 


22.971 


19.309 


16.348 


39.151 


39.003 


37.902 


28.171 




2.0 


22.650 


20.584 


17.311 


14.666 


35.101 


34.970 


33.989 


25.234 






5.0 


20.412 


18.553 


15.621 


13.252 


31.722 


31.607 


30.730 


22.805 




2 


10.0 


19.335 


17.575 


14.808 


12.571 


30.118 


30.010 


29.182 


21.686 




0.0 


41.826 


38.408 


32.776 


27.953 


68.585 


68.377 


66.757 


54.742 






0.5 


28.234 


25.928 


22.160 


18.933 


46.490 


46.357 


45.283 


37.338 




20 


1.0 


24.608 


22.593 


19.324 


16.530 


50.619 


40.506 


39.579 


32.689 




2.0 


21.854 


20.064 


17.190 


14.735 


36.253 


36.157 


35.346 


29.260 






5.0 


19.387 


17.809 


15.318 


13.185 


32.503 


32.427 


31.732 


26.409 






10.0 


18.115 


16.651 


14.370 


12.410 


30.636 


30.573 


29.943 


25.056 



Skew 


Frequency 


angle, ip 




0° 


Mode 1 




Mode 2 


15° 


Mode 1 




Mode 2 


30° 


Mode 1 




Mode 2 


45° 


Mode 1 




Mode 2 



Table 10: Influence of skew angle on linear frequencies for a simply supported FGM plate 
{a/b — l,a/h — 10) with a center crack (c/a ~ 0.2) and with a temperature gradient, 
(Te = 600K, Trn ^ 300K) 

gradient index, n 

0.0 0.5 1.0 2.0 5.0 10.0 

16.605 11.273 9.859 8.799 7.788 7.420 

41.754 28.467 24.944 22.350 20.187 19.150 

18.351 12.480 10.918 9.751 8.748 8.257 

41.778 28.487 24.962 22.365 20.202 19.169 

23.956 16.341 14.310 12.796 11.519 10.911 

47.407 32.365 28.365 25.416 22.969 21.814 

35.793 24.473 21.448 19.200 17.325 16.450 

61.697 42.180 36.978 33.134 29.956 28.477 



The influence of thermal gradient on the linear frequencies of FGMs is 
examined in Tables [8] and [9] and graphically in Figure ([8]) with different 
surface temperatures. The temperature field is assumed to vary only in the 
thickness direction and is determined by Equation Q. The temperature for 
the ceramic surface is varied (T^ = 400K, 600K) while maintaining a constant 
value on the metallic surface (T^ = 300i^). Figures IQ and f lTO|) shows the 
first two mode shapes of simply supported FGM plates with and without 
a crack. A similar trend can be observed for the FGM plates in thermal 
environment. The natural frequency decreases with increasing gradient index 
and crack length. Also, the frequency further decreases with increase in 
temperature gradient as expected. The mode shape (second mode) for a 
simply supported square plate with a center crack is shown in Figure ( ITTj) 
for different crack length and a/h = 10, Tc = 600, T^ = 300, n = 0. It is 
seen that with increase in crack length, the plate becomes locally flexible. 

Lastly, the effect of skewness of a simply supported and clamped square 
FGM plates in thermal environment is investigated. The thickness of the 
plate is taken as a/h=10 and the plate is assumed to have a center crack 
c/a = 0.2. The results are shown in Tables [TO] and [TT] for various gradient 
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(a) Mode 1 




(b) Mode 2 

Figure 9: Mode shape for a square plate without a crack. The mode shape is normahzed 
with the maximum value. 
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(a) Mode 1 




(b) Mode 2 

Figure 10: Mode shape for a simply supported square plate with a crack (cj a — 0.2, ajh = 
10; Tc = 600, T„, = 300, n = 0). The mode shape is normalized with the maximum value. 
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(a) c/a=0.2 




(b) c/a=0.4 




(c) c/a=0.6 



Figure 11: Mode shape (second mode) for a simply supported square plate with different 
crack sizes {a/h = 10, Tc = 600, Tm = 300, n = 0). The mode shape is normahzed with 
the maximum value. qn 



Table 11: Influence of skew angle on linear frequencies for a clamped FGM plate [a/b - 
l,a/h = 10) with a center crack {c/a — 0.2) and with a temperature gradient, {Tc : 
600K, Tm = 300K) 

gradient index, n 

0.0 0.5 1.0 2.0 5.0 10.0 

29.331 19.991 17.501 15.666 14.147 13.431 

57.569 39.137 34.283 30.691 27.726 26.335 



Skew 


Frequency 


angle, ip 




0° 


Mode 1 




Mode 2 


15° 


Mode 1 




Mode 2 


30° 


Mode 1 




Mode 2 


45° 


Mode 1 




Mode 2 



30.963 


21.110 


18.482 


16.546 


14.944 


14.189 


57.338 


39.178 


34.318 


30.724 


27.755 


26.382 


36.649 


25.010 


21.903 


19.611 


17.717 


16.831 


63.181 


43.191 


37.834 


33.865 


30.590 


29.082 



49.730 33.987 29.773 26.658 24.088 22.899 
78.840 53.951 106.69 42.281 38.179 36.311 



indices and skew angles. It can be seen from Tables [TO] and [TT] that with 
increase in skew angle ip, the linear frequency increases and with increase in 
gradient index n, the linear frequency decreases. 

5. Conclusion 

The linear free flexural vibration behavior of FGM plates with and with- 
out thermal environment is numerically studied using QUAD-8 shear flexible 
element. The formulation is based on first order shear deformation the- 
ory. The material is assumed to be temperature dependent and graded in 
the thickness direction according to the power-law distribution in terms of 
volume fractions of the constituents. The effective material properties are es- 
timated using Mori-Tanaka homogenization method. Numerical experiments 
have been conducted to bring out the effect of gradient index, aspect ratio, 
crack length, thickness and boundary condition on the natural frequency of 
the FGM plate. From the detailed parametric study, it can be concluded 
that with increase in both the gradient index and crack length, the natural 
frequency decreases. In both cases, the decrease is due to stiffness degrada- 
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